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ABSTRACT 

Sagittarius A*, the supermassive compact object at the center of the Galaxy, exhibits outbursts in the near 
infrared and X-ray domains. These ^are,? are likely due to energetic events very close to the central object, 
on a scale of a few Schwarzschild radii. Optical interferometry will soon be able to provide astrometry with 
an accuracy of this order (~ 10 /las). In this article we use recent photometric near infrared data observed 
with the adaptive optics system NACO at the Very Large Telescope combined with simulations in order to 
deploy a method to test the nature of the flares and to predict the possible outcome of observations with the 
Very Large Telescope Interferometer. To accomplish this we implement a hot spot model and investigate its 
appearance for a remote observer in terms of light curves and centroid tracks, based on general relativistic ray 
tracing methods. First, we use a simplified model of a small steady source in order to investigate the relativistic 
effects qualitatively. A more realistic scenario is then being developed by fitting our model to existing flare 
data. While indications for the spin of the black hole and multiple images due to lensing effects are marginal 
in the light curves, astrometric measurements offer the possibility to reveal these high-order general relativistic 
effects. This study makes predictions on these astrometric measurements and leads us to the conclusion that 
future infrared interferometers will be able to detect proper motion of hot spots in the vicinity of Sagittarius A*. 

Subject headings: astrometry - black hole physics - Galaxy: center - gravitational lensing - techniques: inter- 
ferometric 



1. INTRODUCTION 

The center of our Galaxy is by far the closest galactic nu- 
cleus to Earth, about 100 times closer than the center of 
Andromeda. The radio source Sagittarius A* (Sgr A*) lies 
motionless (with a motion constraint of less than 45 km/s, 
iReid et al1l2007h at the center of gravity of the nuclear clus- 
ter. It has variable counterpart s in the in frared an d X-ray do- 
mains (Bagano ff et al. 2001; Porqu et et al. 2003; Genzel et alJ 
l2003t lEckart et al. 2004). Dynamical evidence shows that 
this source is powered by a supermassive object, most likely 
a black hole (jS chodel et al. 2002, 1200I iGhezetalJ llooa 
lEisenhauer etani2005ai) . Its mass Mq = 3.61 ±0.32 x IO^'Mq 
and distance Rq = 7.62 ± 0.32 kpc can be deduced from the 
motion of the in nermost stars, which are d ominated by a point 
mass potential (lEisenhauer et alJ l2005al) . The correspond- 
ing Schwarzschild radius is given by rs = 2QM/c^ ~ 15 Rq, 
yielding an apparent size of 6^ = rs/Ro ~ 9.3 /ias on sky. 
In contrast, for a stellar black hole (M ~ IM0) at 1 kpc, 
9s —2x 10"^ /ias. Indeed, Sgr A* is the black hole candidate 
with the largest apparent size, and therefore the candidate of 
choice for resolving effects at the Schwarzschild scale. 

In the near infrared (NIR), Sgr A* is studied using 10m- 
class telescopes such as the ESO Very Large Telescope (VLT) 
and the Keck telescopes with a diffraction limit of a few 
10 mas. Both the VLT and Keck feature several telescopes 
on the same site (4 and 2, respectively). They offer an inter- 
ferometric mode in which the light collected by two or more 
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of these telescopes can be combined in a coherent fashion, 
allowing phase-referenced measurements to be conducted. 
Used with a proper data acquisition and reduction schemes, 
these interferometers can deliver reconstructed images with 
resolutions as good as a few mas. Even in the case of Sgr A*, 
this is still two orders of magnitude larger than rs. However, 
these interferometers can also provide differential astrometry 
between the target and a reference source with a very high ac- 
curacy. The Phase-Referenced Imaging and Micro-arcsecond 
Astrometry (PRIMA, Glindemann & Leveque 2000) facility 
on the VLT Interferometer (VLTI) will reach 10 /ias accuracy 
for a typical observational setup: f = 30 min integration time, 
a separation of = 10" between the reference star and the 
science object, and a (projected) baseline b = 200 m. This 
accuracy scales as dShao & Colavita 1992j): 
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We have started studying the possible outcome of high- 
accuracy astrometry of matter very close to Sgr A* as 
part of the PRIMA R eference Missi ons ESO documen t 
(iDelplancke et al.ll2004 . Further study (iPaumard et al.ll2005h 
has led us to propose a second-generation instrument for the 
VLTI: GRAVITY, a general purpose instru ment optimized 
for the Galactic Center (GC) science cases (Eisen hauer et alJ 
2005b). The GRAVITY-consortium has completed a phase A 
study for this instrument, which was recommended to ESO by 
the scientific advisory committee in late 2007. 

1.1. Observational facts about Sgr A * 

Sgr A* is v ery faint compared to its Eddington luminosity 
(L < 1 0~^Lfhh ■ iBaganoff et al.]l2003h . However, it shows vari- 
ability at various timescales. In particular, it exhibits outbursts 
of en ergy called flares in both X-rays (e.g. Baganoff et aD 
2001L l2003h and NIR (e.g. iGenzel et al.1 l2003flGhez etall 
2004l:IClenet et al.ll2005h . These flares typically last for a few 
hours and occur several times a day. On these occasions, the 
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luminosity of Sgr A* increases by a factor of up to 100 in X- 
rays and up to 20 in K-band, making Sgr A* a niK ~ 15 source 
jGenzel et alj |2003|) . Superposed upon these bright, hour- 
scale flares, quasi-periodic variations with a typical timescale 
around 20 rnin are observed in the NIR (iGenzel et alJ 120031: 

lEckart e t alJl2006HTrippe et al. 2007). 

Dodds-Eden et al. (in preparation) and lPorquet et al.l ( l2008h 
observed a bright flare from the GC simultaneously in NIR 
and X-ray on April 4, 2007. The power spectrum of the 
NIR light curve gives indication for the existence of a quasi- 
periodic signal at about 23 min. This characteristic timescale 
can be interpreted as the orbital period 




of matter on a circular orbit of radius r around a rotating black 
hole (Kerr black hole with gravitational radius r „ := rs/2) with 
dimen sionless spin-parameter a (a G [-1; 1]) dBardeen et alj 
[1971) . It is a natural timescale, if the flares involve mate- 
rial on these orbits close to Sgr A*. However, this orbit is 
only stable for radii larger than the innermost stable circular 
orbit (ISCO), which itself depends on a (see iBardeen et alJ 
Il972l) . Thus, there exists a lower boundary on the orbital 
timescale for a given spin-parameter In the case of Sgr A*, 
this timescale lies somewhere between 27.3 min (for a = 0) 
and 3.7 min (for a = 1, counter-rot ating orbits with a < are 
not considered here). For instance. iGenzel et alj ( l2003b find a 
quasi-periodic signal at roughly 17 min. Assuming this varia- 
tion corresponds to the orbital period on the co-rotating ISCO 
around Sgr A*, they give a lowe r limit for its spin-parameter 
(a > 0.52). iTrippeetalJ (l2007h argue that flares happen at 
various radii larger than the ISCO. From the shortest period 
ever observed (13 ± 2 min), they argue thata> 0.70. 

In addition to this, lEisenhauer et al.l (12005a ') were able for 
the first time to derive the slope of the spectral energy distri- 
bution (SED) of GC- flares in the K-band. T hey found a red 
spectrum: i^Si, oc v'^ . iGillessen et al.l (l2006l) showed in addi- 
tion that the spectral index varies significantly during a flare 
and is correlated with the concurrent brightness. While dim 
states show a red spectrum, the fla re appears incre asingly blue 
as its brightness rises. In contrast. iHornstein et al.l (120071) find 
a constant SED-slope for the flaring states of Sgr A* with a 
spectral index being close to zero. In any case, these observa- 
tions indicate the flares being due to synchrotron emission. 

Recent observations of the GC in the NIR also included 
polarimetric measurements of Sgr A * dEckart et al.l 120061: 
iTrippe et ani2007h . iTrippe et aP (l2007h find significant hnear 
polarization in the flares with polarization degrees of 15%- 
40%. While they see a correlation between the variability of 
polarized flux and total flux, the degree of polarization seems 
to remain constant on the orbital timescale. Furthermore, an 
increase of the polarization degree and a swing of the po- 
larization angle of roughly 70° within 15 min is reported in 
one case. It happens in the declining phase of the flare. This 
could be interpreted in a scenario where the emitting material 
leaves the plane in which is was orbiting, perhaps through a 
jet. These observations strongly favor the nature of the emit- 
ted light as synchrotron radiatio n and that the accretion d isk 
is seen nearly edge-on (compare iBroderick & Loebll200 6b'). 

In the radio domain, Sgr A* appears as a permanent, 
resolved source. Scattering effects in the interstellar dust 
ma ke it appear larger, w ith an apparent size pro portional to 
dPavies et all 1 19761) . iBower et all d2004 and IShen et all 



d2005h were able to measure the intrinsic size of the emission 
region at 3.5 and 7 mm, finding 0.126 ±0.017 mas ~ 15rs and 
0.268 ± 0.025 mas ~ 30rs, respectively (maj or axis). The lat- 
est measurements were carried out by Doeler nan et alj d2008h 
at 1.3 mm that set a size of about 37 /ias on the intrinsic di- 
ameter of Sgr A*. Hour-s cale variability has been ob served 
at 2 and 3 mm jMiyazaki et al.ll2004l: iMauerhan et"an r2005). 
lYusef-Zadeh et al.. 42006) have detected flaring activity at 
7 mm and 1 3 mm, with a 20^0 min delay between the peaks 
of the two lightcurves. At centimeter wave lengths, periods 
of 57 and 106 days have been claimed (Townes et al.lll999t 
Zhaoet al.H 2001). but were later refuted ( Macquart & Bowej 
20061) . Polarimetric measurements yielded polarization frac- 
tions around 10% and hig hly variable polarization angles on 
timescales of a few hours dBower et alJl2005al : lMarrone et alj 
l2006l) . 

At the currently possible imaging resolution (~ 1" for 
Chandra, 60 mas for the VLT in the K-band), the flares oc- 
cur at the same location as the unresolved X-ray quiescent 
emission. This position is consistent (within the error bars of 
a few mas or a few hundred rs) with the focus o f the orbits 
of nea rby stars orbiting the black hole (Schodel et alj|2002l 
I2003t iGhez et all l2003t lEisenhauer etali ,2005a). It is also 
coincident with t he radio source, which appears motionless 
dReid et al.ll2007 l). However, on shorter timescales (between 
50 and 200 min) a centr oid wander of about 100 /ias has been 
observed by iReid et al.l ((2008) with the Very Long Baseline 
Array (VLBA). From their observations they can rule out hot 
spots with orbital radii larger than 7.5rs. 

The purpose of this paper is to demonstrate the feasibility 
of measuring the position and proper motion of NIR flares 
with an accuracy of 10 /ias. In order to resolve this proper 
motion, one must not only reach the necessary astrometric 
accuracy, but also resolve flares temporally. Because the sam- 
pling must exceed the characteristic periodicity seen in the 
NIR light curves, this represents a very challenging experi- 
ment. 

In the typical observational setup mentioned above, this ac- 
curacy requires a 30 min integration time in order for atmo- 
spheric residuals to be averaged out. Given Eq. [T] the re- 
quired integration time for a fixed accuracy depends on the 
separation Q. In the GC a suitable niK = 10 reference star 
can be found as clo se to Sgr A* as 1.21" (GCIRS 16NW, 
iPaumard et alJl2006h . With such a close reference object and 
a projected baseline of 100 m, the required integration time 
for 10 /ias astrometry reduces to roughly 1 min. Therefore, 
the VLTI infrastructure is able to provide time-resolved as- 
trometry of GC-flares at a sufficient accuracy. However, this 
is a very demanding experiment and requires a specifically 
optimized instrument. The GRAVITY experiment has been 
propose d to fulfill all the requirements for the GC science case 
(Eisenhauer et al. 2005b). 

Similar attempts for the astrometric o bservation of S gr A* 
are being planned for the radio domain dBoweiil2005bl) . The 
VLBA is expected to improve its imaging resolution down 
to 20 /ias within the next decade (Doeleman & Bower 2004). 
However, for technical reasons the integration times in ra- 
dio and submillimeter wavelengths are significantly longer 
than in the NIR. Thus, no dynamic picture will be accessi- 
ble easily, rather a time-averaged image can be generated. 
This has already been stu died by means of simulations by 
IBroderick & Loebl d2006b). Generally, enhanced opacity at 
longer wavelengths smoothes out the strong features in the 
light curves and centroid tracks expected in the NIR. 
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1.2. Flare-models 

Several hypotheses have been proposed to explain the ori- 
gin of the flares from the GC. The most common ones can be 
summarized as follows: 

A Star-disk interac tions rather far away from the black 
hole (10^-lOVQ. lNavakshin et al]l20M) : 

B Synchrotron emission from electrons accelerated in the 
central part (< 10 rs), through processes such as turbu- 
lent acceleration , reconnections, and weak shocks (e.g. 
lYuan etani2004h : 

C Sudden heating of hot electro ns (e.g. by magnetic re- 
connection) in a permanent jet dMarkoff et alj|20()lh : 

D Increase in the accretion rate due to the infall of a 
clump of material, being rapidly driven inwards due to a 
Rossby wave instability (RWI) dTagger & Meliall2006l 
iFalanga et al...2007 ). 

The star-disk interaction model appears unlikely. It is un- 
able to explain either the observed synchrotron emission or 
the quasi-periodic signal on a timescale of minutes. More- 
over, actual constraints on Sgr A*'s disk size in the Radio 
(/"disk ^ 1 At/. IShenl i2008) are in conflict with the assumption 
of stars remaining stable on such close orbits. 

No matter what the emission process producing the flares 
is, the matter responsible for it must move with at least the 
order of the Keplerian circular velocity 

[gM~o c [t; ( rY'^ ,^ . _i 

= —=\ — ~ — xlOuasmm (3) 
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This is true irrespective of whether this matter is orbiting the 
black hole, falling into it in free fall, or escaping in a jet or 
wind. Given that r < 1 mas ~ lOOrs, we have v > 1 /lasmin"'. 
Therefore, this matter moves by several times the proposed 
accuracy (10 /ias) during one flare. This typical velocity is 
large compared to the apparent velocity of the star proposed 
as interferometric phase reference. 

If the flares are due to magnetic reconnection in a jet (sce- 
nario C), sufficiently high above the accretion disk, the elec- 
trons responsible for the emission would have relativistic ve- 
locities, but the emission region may not move at the same 
speed. This region is tied to the node in the magnetic field 
where the reconnection happens. If this node drifts on the 
hour timescale, some long term proper motion could be ob- 
served. In this scenario, the short-term variability in the 
lightcurves could come from individual blobs of denser ma- 
terial being ejected at a pseudo-frequency, typical for the dy- 
namical timescale of the inner disk. When passing through 
the node at a significant fraction of c ~ 15 /ias min~', a short 
term proper motion could be expected from these individual 
blobs. 

Ho wever, the observational data at hand (see iTrippe et al.l 
l2007l for a summary) strongly favor another interpreta- 
tion, namely the emission of synchrotron radiation from 
matter orbiting the black hole near its ISCO. This matter 
would be confined to a small region in the disk (r < 0.3rs, 
iGillessen et al.ll2006l) containing highly accelerated electrons 
in a broken pow er-law distribution {n^ w 10^' cm"^7"^ in the 
NIR, lYuan et al. ''20b4). Such a radiating blob could be pro- 
duced by scenario B. Despite the very high velocity of the 
material (about half the speed of light), no long-term proper 



motion would be detected at the scale of more than one orbital 
period. Nevertheless, this scenario opens up extremely inter- 
esting scientific perspectives. Time-resolved astrometry of the 
flares would allow the shape of the ISCO to be traced, lensed 
by the strong gravitational field of the black hole. This would 
give access to the orientation of the orbit as well as other phys- 
ical parameters of the model, and more fundamentally provide 
a probe for the strong gravitational field of a black hole. 

This scenario is akin to the RWI scenario (scenario D). 
In both cases, the emission originates ve ry close to the 
black hole, within a few Schwarzschild radii ( Tagger & Melial 
l2006l) . However, there is a fundamental difference: instead of 
a blob, sc enario D produces a somewhat symmetrical spiral 
structure jFalanga et al.ll2007h . It is created by a RWI in the 
accreting plasma producing quasi-periodic oscillations in the 
accretion rate which are close to the orbital period. Though, 
these oscillations are not of Keplerian origin and therefore 
do not allow the black hole spin to be constrained from the 
periodic timescale. They are caused by a pattern rotation in 
the disk, a wave phenomenon. Further studies are needed to 
assess whether this process could mimic the observed light 
curves and whether it could lead to an observable astrometric 
signature. 

2. SIMULATIONS 

We will hereafter always assume that flares from Sgr A* are 
due to matter on Keplerian orbits around the black hole and 
that the periods observed in the light curves are a signature of 
these orbits. After a discussion of our basic methods we will 
first investigate the consequences of this scenario on the basis 
of a simplified model. Therefore, we study the appearance of 
emission from a small sphere on the co-rotating ISCO about 
a Kerr black hole and generate photometric, astrometric as 
well as polarimetric data. We will then address an advanced 
model with an extended emission region as a phenomenolog- 
ical model for a variable source. 

The lensing effects in the strong gravitational field of a 
black hole and the optical appearance of accretion disks and 
close-by objects ha ve been computed a number of times 
(e g.p unningham & Bardeenlll 9731: iLuminetlll 9791: iBao et aP 
[l994rTKUng et al.i |200(?). Derivations of analytical lens- 
equations, giving both the position of multiple images of a 
source and the time of arrival of each of these images have 
also been attempted. These purely analytical studies are how- 
ever limited in the computation of higher order images (no 
more than second order) and do not account for extended 
sources with variable intrinsic brightness. 

The apparent path on sky of a particle revolving around a 
black hole in a circular orbit of radius r {r = 3rs for the ISCO 
of a Schwarzschild black hole) does not have the shape of 
an ellipse due to strong deflection of light rays by the grav- 
itational field. Observed in high inclination, the orbit's near 
side appears as half-ellipse (the usual appearance of an in- 
clined circle), whereas the far side is de-projected and closely 
resembles a semi-circle. 

In addition, light rays emitted by a nearby source are bent so 
strongly that this source appears at several different locations 
on sky. In principle, infinitely many such images emerge. 
However, higher order images are exponentially suppressed 
in angular diameter and have increasingly longer light travel 
times compared to the primary image (see Fig. IB- Moreover, 
absorption increasingly diminishes the brightness of images 
with longer light paths. 
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2.1. Phenomenological assumptions 

Unlike the Newtonian case, in both Schwarzschild and Kerr 
geometry bound orbits are not closed in general. Eccentric 
orbits in the Schwarzschild metric or in the equat orial plane 
of the Kerr metric are planar, but not closed (see iBao et all 
|1994|) . Due to Lense-Thirring precession, circular orbits out- 
side this eq uatorial plane are not planar and can be described 
as helices ('Wil tdnslll972[) . For this reason, a cloud falling 
onto the black hole on a trajectory inclined with respect to the 
equatorial plane would soon form a geometrically thick disk, 
symmetrical with respect to the equatorial plane, and settle on 
this plane to be circularized. In the following we will there- 
fore concentrate on circular orbits in the equatorial plane. 

The power law SED observed in the GC-flares with 
a red spectral index i n the NIR ([Eisenhauer et"an l2005al : 
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iGillessen et al] 120061; iKrabbe et al. 
chrotron emission. This indicates the phenomenon to be dic- 
tated by a strong magnetic field. Po ssible c onfigurations of 
such magnetic fields are discussed by lUzde nsky (2005). The 
field lines should intersect the accretion disk orthogonally at 
its surface. We will therefore consider the field to be poloidal 
near the disk. On the other hand, the magnetic field inside the 
accretion disk is likely frozen. The field lines must warp and 
align with the velocity field, leading to a toroidal geometry in- 
side. Shear tends to make the magnetic field lines loop around 
the black hole, leading to magnetic reconnection events. The 
stronger shear near the inner edge of the disk causes more 
reconnection events towards this edge. We will consider a 
bright sphere close to the ISCO, which should be a good ap- 
proximation for the inner edge of the quasi-circular accretion 
disk. The sphere is just a first approximation of the emitting 
region geometry. Its real shape presumably follows the field 
lines and would therefore extend as an arc along the orbit. 

2.2. Relativistic my tracing 

We choose to use relativistic ray trac i ng to model t he light 
bending effects (IWeiskopfll2000l l200lt lMulledl2006h . Light 
rays are represented by null geodesies which are determined 
by integrating the equation of motion for photons in Kerr 
spacetime. Initial values for this ordinary differential equa- 
tion are determined by the position, orientation, and field of 
view of the observer's camera and by the coordinates of the 
corresponding pixel in the image plane. Starting from the ob- 
server, the geodesies have to be integrated back in time until 
they fulfill one of the breaking conditions: they either reach 
the surface of the source (as in Fig. [T]!, they cross the hori- 
zon of the black hole, or they leave a certain bounding box 
around the region of consideration. A further limit is given by 
a maximum length for any computed geodesic. 

The main advantage of this method is that it does not yield 
any analytical approximation, it is exact up to the numeri- 
cal precision. The fact that different null geodesies do not 
affect each other allows the code to be parallelized, i.e. the 
field of view can be divided into smaller sections to be com- 
puted on several CPUs and reassembled later on. We set up 
a very simplified configuration of a nonrotating sphere of ra- 
dius 0.25rs on a circular orbit in the equatorial plane of a Kerr 
black hole. The sphere's surface is homogeneous and emits 
light isotropic ally. From the ray tracing code we compute the 
lensed images of the source (we typically resolve images up 
to the fourth order) with the corresponding light travel times 
t and frequency shifts z := t'obs/^'em computed for each pixel 
('obs' for observed and 'em' for emitted). In addition, the 




Fig. 1. — Travelled light paths of the primary (solid), secondary (long- 
dashed), tertiary (short-dashed) and quaternary (dotted) image of the source 
at a given anival time (t = 0) at the observer's location. The higher the order 
n of an image I„ , the earlier its emission time ; and the smaller the absolute 
value of its deflection angle 

code provides us with the local tetrad Cq of the observer, par- 
allel transported to the source in each pixel. It is necessary for 
the polarimetric analysis of the flare model. 

The main restriction to the ray tracing method is a tradeoff 
between resolution and computing time. In order to keep the 
computing time manageable, we had to consider frames with 
no more than 1000 x 1000 pixels. In addition to this, again for 
the sake of reasonable computing time, it has been necessary 
to use a finite distance between the observer and the black hole 
(60rs). However, this yields a good enough approximation, 
since we are able to resolve all multiple images up to the third 
order. 

For our simulations we set up a grid of four inclinations 
(/ = 20°, 50°, 70°, 90°), four orbital radii (r/nsco = 1.0, 1.2, 
1.5, 2.0) and four black hole spin-parameters (a = 0, 0.52, 0.7, 
0.998). All possible combinations yield a grid of 64 different 
setups, allowing to analyze the impact of these parameters. 
For each configuration we ray traced one complete orbit, con- 
sisting of roughly 100-500 frames depending on the orbital 
period of the source (for the a = 0.998 case only two simula- 
tions with r = 2.0risco were carried out). 

2.3. Light curves and centroid tracks 

Assuming the source to emit a power law SED in its rest 
frame in the NIR, 

z/e„,5r (X i^:'^ (4) 

with spectral index a, we can compute the apparent surface 
brightness 5°''^ (as seen by a distant observer) from the fre- 
quency maps provided by the ray t racing code. We simply 
have to apply Liouville's theorem (IMisner et al.lll973b . stat- 
ing: 

S7/i^L = Sfy>yL = const. (5) 
For the apparent surface brightness this yields 

cobs 4 _ . 



^obs 



(6) 



depending on the generalized redshift factor z, which is re- 
sponsible for relativistic beaming- and Doppler effects, as 
well as gravitational redshift. The result is that the source 
appears brighter and blueshifted on the approaching side, re- 
spectively fainter and redshifted on the receding side. We are 
not interested in the absolute value of since we do not 
model the emission process in a self-consistent manner We 
only study the dependencies of the simulated light curves on 
the various parameters of the model. 

Light curves are created by integrating the apparent sur- 
face brightness over a certain frequency range and over the 
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field of view: 

L:=yy 5f di^obsdAfov (7) 

Thus, every simulated frame produces one data point of the 
light curve. We normalize each light curve by its minimum 
flux and plot its magnitude (m = -2.5121og[QL) against the 
orbital phase (p in radians (0 = designates the point where 
the hot spot appears closest to the observer). 
Similarly, the centroid of emission is computed as 

C:=iyy"r5f dz/obsdAfo, (8) 

with r being the vector from the origin (center of field of view) 
to the pixel position. The centroid track is generated by plot- 
ting C for every frame. We always use the mass and distance 
of Sgr A* for the scaling of its apparent size (Irs corresponds 
to about 9.3 /ias). 

3. RESULTS 

3.1. Small, constant source 

First we study the observable properties of a constantly ra- 
diating sphere of radius 0.25rs and spectral index a = on 
a co-rotating circular orbit around a generally rotating black 
hole. The sense of rotation of the black hole (and thus the 
direction of rotation of the source) is arbitrarily chosen to be 
counterclockwise. This approach is suited to investigate the 
main features of the hot spot model. 

3.1.1. Photometry and astrometry 

Orbital radius and inclination — Fig. |2] displays light curves 
and centroid tracks of this model for different orbital radii and 
inclinations of the hot spot. Decreasing the orbital radius gen- 
erally leads to narrower, higher peaks in the light curves. This 
is due to the increasing spot- velocity leading to stronger rela- 
tivistic Doppler- and beaming effects. The same changes can 
be produced by increasing the inclination of the orbit, since 
this increases the line of sight velocity of the hot spot. In 
our chosen parameter range, changes in the inclination ex- 
hibit these relativistic effects more strongly than changes in 
the orbital radius. 

The light curves in 90° inclination reveal an additional fea- 
ture: before each of them reaches its usual peak-maximum 
(approaching node), a spiky second peak appears. Its peak- 
flux exceeds the usual maximum by far It originates from 
the development of an Einstein ring, an apparent ring-shaped 
image of the hot spot with the black hole in its center It oc- 
curs when the light source, the black hole and the observer are 
aligned (and the black hole is in between). Due to the finite 
extension of the source, conditions for an Einstein ring are ful- 
filled in a finite (although small) range of inclinations. Thus, 
it is no mathematical artifact of the exact 90°-case. Higher 
order Einstein rings appear, when light rays revolve the black 
hole. For instance, a secondary Einstein ring appears when 
the hot spot passes in front of the black hole {(f) = 0). We then 
receive photons from the far side of the spot, which have or- 
bited the black hole once. In fact it occurs a bit later, since one 
has to consider the finite speed of the photons when orbiting 
the black hole. The corresponding peaks can also be seen in 
the 90°-plot (at (f) « 7r/2), although they are a lot fainter than 
the ones caused by the primary Einstein ring. 

At lower inclinations we no longer expect Einstein rings, 
but still multiple images of the hot spot. Obviously, their ef- 
fects on the light curves are marginal. Only the light curves 



with 70° inclination show a slight distortion of the main peak, 
suggesting the existence of a secondary image. It is not only 
its small apparent size that prevents the secondary image from 
being observable in the light curves, but also the fact that its 
brightness is not as variable as that of the Einstein ring. 

Turning to astrometry, decreasing the orbital radius obvi- 
ously leads to smaller centroid tracks. Then, as with the light 
curves, the relativistic effects increase, i.e. the deviation of the 
centroid tracks from the Newtonian case (projected circles) 
becomes more significant. On the one hand, this deviation 
is caused by multiple images of the hot spot. When the sec- 
ondary image gets brightest, the centroid of the received light 
lies somewhere in between the primary and the secondary im- 
age. But since the latter only brightens for a short time (it 
moves faster on the approaching side), this causes a cusp in 
the centroid track. On the other hand, deflections of the pri- 
mary rays result in a further distortion of the centroid track. 
The rear part of the orbit appears to be folded upwards at 
higher inclinations. 

Again, raising the inclination leads to stronger relativistic 
effects. The cusp in the centroid track, for instance, becomes 
more and more prominent. As shown in the lower right plot of 
Fig.|2] the classical sine-pattern of the 90°-case is deformed. 
It is stretched when the source is receding and compressed 
when approaching (due to the finite propagation speed of the 
emitted light rays). In addition, distinct spikes appear in the 
pattern. These can again be associated with the Einstein rings. 
Comparing the occurrence of these spikes with the occurrence 
of the sub-peaks in the light curves, one can ascertain that they 
appear at the same orbital phase. 

Spectral index — Changing the SED of the hot spot emission 
has strong implications for the light curves as well. Since we 
are observing in a confined range of frequencies (observing 
band), Doppler shifting the emitted radiation makes us ob- 
serve different parts of the SED every instant of time. If, for 
instance, the hot spot is approaching us, the observed radi- 
ation is blue-shifted. Depending upon whether the hot spot 
emission follows a blue, white or red SED, the observed flux 
is lower, equal or higher, respectively. This means that a red 
SED {a < 0) amplifies the beaming effect by brightening the 
hot spot on the approaching side and dimming it on the re- 
ceding side, while a source with a blue SED (a > 0) behaves 
contrariwise. For a white SED (a = 0) the Doppler shift does 
not change the observed flux at all. 

Fig.[3]displays light curves and centroid tracks for all three 
cases. Clearly, the light curve peaks rise with decreasing spec- 
tral index. In the case of a blue SED the observed flux of 
the hot spot at the approaching side is diminished to such an 
extent that the secondary image of the blob is revealed. Its 
brightness exceeds the total flux of the primary image for a 
short fraction of the orbital period. In the case of a red SED 
the opposite happens: the primary peak gets so narrow and 
bright that the secondary peak appears slightly before the hot 
spot reaches the approaching node. 

For the centroid tracks, it appears that those with a red spec- 
trum are dragged towards the approaching side, whereas those 
with a blue spectrum are more centered towards the reced- 
ing side. The multiple images have a greater influence in the 
case of a red SED, since the primary image of the hot spot 
is dimmed more strongly on the receding side of the orbit, 
whereas the secondary image is amplified more strongly on 
the approaching side. This even makes the cusp from the sec- 
ondary image turn into a loop; the tertiary image may create 
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Fig. 2. — Light curves (top) and centroid tracks (bottom) of a compact hot spot orbiting a Schwarzschild black hole in four different orbital radii and inclinations 
(f/i'isco = 2.0 (solid), 1.5 (long-dashed), 1.2 (short-dashed), 1.0 (dotted), inclinations are given in the upper left corner of each plot). Since the ordinary centroid 
track for the 90° -case just yields a horizontal line, we instead plot the orbital phase </<(0 against x(t). 
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Fig. 3. — Light curves (left) and centroid tracks (right) of a compact hot 
spot with different spectral indices {a = -3 (dotted), a = (solid), a = 3 
(dashed)) orbiting a Schwarzschild black hole in 70° inclination on the ISCO. 

another cusp. The centroid tracks of hot spots with a blue 
SED are generally more symmetric and show less relativistic 
deflection. 

Black hole spin — The impact of different black hole spin- 
parameters on light curves and centroid tracks is shown in 
Fig.|4] Increasing the spin-parameter essentially produces the 
same changes as decreasing the orbital radius of the hot spot, 
both in the light curves and the centroid tracks. This is due to 
the fact that we specify the orbital radius in units of the ISCO, 
which shrinks with higher black hole spin (risco ~ 3rs, 2. Irs, 
1 .7rs, 0.6rs for a = 0, 0.52, 0.7, 0.998, respectively). By look- 
ing more closely at the Einstein ring peaks in the 90° -light 
curves, one can find a slight broadening with higher spin. This 
is because the source occupies a bigger fraction of the orbit at 
smaller radii. The result is that the required conditions for an 
Einstein ring are met for a larger orbital phase (this can al- 
ready be seen in Fig. |2]l. In the high spin case (a = 0.998) the 
shape of the light curves is dominated by multiple images. 

With higher spin-parameters the centroid tracks show even 
greater deformation, i.e. the cusp is enhanced and the rear part 
of the track is dragged up further. In addition, a second cusp 
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Fig. 4. — Light curves (top) and centroid tracks (bottom) of a com- 
pact hot spot orbiting a Kerr black hole with four different spin-parameters 
(a = (solid), 0.52 (long-dashed), 0.7 (short-dashed), 0.998 (dotted), incli- 
nations are given in the upper left comer of each plot, the orbital radius is 
r = 2.0risco)- 

emerges on the receding side, due to the tertiary image of 
the spot. The cusps may again turn into loops provided the 
multiple images are bright enough. The lower right plot in 
Fig.Elshows the 90°-case. Just like the sub-peaks in the light 
curves, the spikes in the deformed sine-pattern are broadened 
with higher spin-parameters of the black hole. In the high 
spin case, Einstein rings dominate throughout the entire or- 
bital motion. Thereby the motion of the hot spot is almost 
completely shrouded. 



7 







□ 


1 1 1 1 1 1 1 1 1 1 1 




[6dlu] uoi 




□ 

A 


D 


□ 


o 
o 

:^ 4 

Q_ 

E 

O 

a; 
> 

Z) 

o 

■1 2 






^ OO 


A 


n 






1 1 1 1 1 1 1 1 1 


, 1 


20 


40 50 80 100 
Centrold deflection [/j-qs] 


120 



Fig. 5. — Amplification of the liglit curve plotted against the horizontal 
diameter of the centroid track for each simulated grid-configuration of the 
compact hot spot. The symbol-shapes indicate different inclinations (; = 20° 
(star), 50° (diamond), 70° (triangle), 90° (square)). The symbol-size indi- 
cates the orbital radius (r/risco = 2.0 (large), 1.5 (medium), 1.2 (small), 1.0 
(tiny)). The symbol-thickness indicates the spin-parameter (ci = (bold), 0.52 
(thick), 0.7 (medium), 0.998 (thin)). 

Aside from its influence on the radius of the ISCO, the black 
hole spin has only marginal impact o n either the shape of light 
curves and centroid tracks (see also iBroderick & Loebll2005l 
l2006ah . For this reason it is presumably impossible to de- 
duce the spin-parameter from the light curves alone. How- 
ever, including polarimetric properties of this model allows 
making stronger statements about the spi n parameter. On th e 
basis of observation al data presented bv lEckart et alj (l2006l) . 
iMeveret alj (l2006allb ) were able to constrain the spin param- 
eter of Sgr A* to the range 0.4 < a < 1 on a 3cr level. 

Another possible approach to measure the spin of Sgr A* 
would be the combination of photometric and astrometric ob- 
servations. This is illustrated in Fig.|5] where we compare the 
maximum-to-minimum flux-amplification in the light curves 
with the maximum horizontal deflection of the centroid tracks 
for each grid-configuration. Obviously, in such a representa- 
tion it is possible to distinguish every single configuration by 
the set of parameters {r,i,a,P}. This is because the incli- 
nation mainly influences the light curve peaks and does not 
change the horizontal centroid deflection that much, whereas 
the orbital radius has a stronger influence on the centroid de- 
flection than on the light curves. Changing the spin-parameter 
of the black hole has no big effect on the location of a config- 
uration with given orbital radius in this plot. Configurations 
with identical locations (i.e. same orbital radius and inclina- 
tion) can be distinguished by their orbital period P, which is 
known from the light curves. With this the spin-parameter is 
determined via Eq. (|2]i. 

3.1.2. Polarimetry 

Recent data sho\y that the NIR flare emission is polari zed 
(lEckart et al.ll2006HMeyer etani2006al: iTrippe et aUboOTl) . It 
is not yet certain whether the magnetic field in the accretion 
disk is mostly poloi dal or toroidal, but there is evidence both 
fro m simulations bylDe Villiers et al.l (l2003h and observations 
by'Trippe et al. (2007) in support of a toroidal configuration. 
Assuming the magnetic field within the hot spot (compact 
sphere) to follow this geometry one can investigate the polari- 



metric properties (temporal evolution of polarization fraction 
and angle) of the rotating hot spot scenario^ 

A detailed treatment was carried out bv IBroderick & Loebl 
(|2006a). They find that the polarized flux follows the evolu- 
tion of the total flux with a short delay. Additionally, small 
dips in the polarized flux light curve are produced by the mul- 
tiple images of the source. The variation of the polarization 
angle is mainly influenced by the inclination of the hot spot 
orbit. Within one orbital period it rotates continuously, but is 
punctuated by rapid rotations (around 90°) at times where the 
dips in the polarized light curve appear In case of an edge-on 
view, the angle remains mostly constant on account of ax- 
ial symmetry, but the polarization fraction drops significantly 
du ring Einstein ring e vents . 

Mever et al.l (l2007l) used polarimetric data together with hot 
spot simulations to estimate the orientation of the Sgr A* sys- 
tem. They investigated two different magnetic field geome- 
tries (sunspot-like and toroidal) and derive the position angle 
of the Sgr A* equatorial plane normal to be between 60° and 
108° (east of north) on a 3(t level. Besides this, a high incli- 
nation angle of the hot spot orbit is favoured. 

In the context of the VLTI we also want to study the effect 
of polarization on the centroid motion of the flares. Astro- 
metric measurements using polarizing filters may help filter 
out physically different components of the emission region 
(e.g. sphere and arc, see section 13.2b and thus increase the 
apparent proper motion of the single components. Further- 
more, the signatures of multiple images may be easier to de- 
tect, since they generally possess different polarization angles 
with respect to the primary image. To this end we have sim- 
ulated centroid tracks produced by a small, constant source, 
observed through two perpendicular polarizing filters in the 
NIR. 

The results are shown in Fig.|6] The effect on the astrometry 
can be very important when the polarization fraction p of the 
source is high. The centroid tracks show pronounced features 
caused by the multiple images of the hot spot, especially when 
observed in the polarization component along the orbital line 
of nodes. However, most of the flux is blocked by the filter 
in this observing mode, as can be seen in the light curves. It 
only rises up shortly when the multiple images (respectively 
Einstein rings) are amplified most strongly, causing the swing 
in the polarizing angle. Thus, an auspicious observation of 
this kind might be hard to achieve. Moreover, observational 
data support a rather low polarization fraction (Tripp e et al.l 
(2007) mostly see 10%-15%). As apparent from Fig. |6l a 
polarization fraction of 50% already reduces the polarimetric 
deviations from the common centroid tracks to less than the 
projected accuracy of order 10 ^as. 

3.2. Extended, variable source 

The conception of a single compact object revolving on a 
close orbit about a supermassive black hole, maintaining a 
constant shape and brightness is somewhat artificial and has 
presumably little physical relevance. In order to improve this 
model we implement, in addition to the emitting sphere, an 
elongated arc that spreads from this sph ere along the orb it. A 
similar model was already adopted by Meyer et al.l (12006 a,b) . 
it is mainly motivated by the broad shape of the observed light 
curve peaks from Sgr A*. By fitting their model to observed 
flare data they find a constant elongation of the arc trough- 
out the flare. This can also be inferred from the fact that the 
width of the major peaks in the flare light curves does not 
change significantly during one flare. A steady broadening of 
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Fig. 6. — Light curves (top) and corresponding centroid tracks (bottom) in polarized liglit, produced by a small sphere on the ISCO around a Schwarzschild 
black hole seen in 45° and 90° inclination. For the source we assume a spectral index of a = 0, a toroidal magnetic field, and a polarization fraction of 100%, 
respectively 50% (polarization fraction and inclination are given in the upper left corner of each plot). The polarizing filter axis is oriented along the orbital line 
of nodes (dotted lines) or along the orbital spin axis (dashed lines). The solid curves indicate the total flux. 



the peaks, for instance, would indicate gravitational shear of 
the emitting region. However, this does not seem to happen 
given the data at hand. 

The crucial assumption for considering gravitational shear 
is to neglect any bonding forces within the plasma around 
the black hole. However, a strong magnetic field is likely 
to be present and may prevent hot spots from being sheared 
on the Keplerian timescale. The characteristic shear-time 
is then increasing with higher magnetic fi eld strengths, due 
to a m agneto-rotational instability (MRI) (iBalbus & Hawleyl 
1199 Ih . It acts against the tidal forces caused by the gravi- 
tational field and creates an attractive force among adjacent 
elements of the accretion disk, trying to enforce a rigid rota- 
tion. Thus, in the case of a very strong magnetic field, the 
hot spot can be regarded as a rigid source, retaining its shape 
during the whole duration of a flare. 

We vary the intrinsic brightness of the two components 
(sphere and arc) with time to emulate heating and cooling 
sequences. We assume that the flares are due to a heating 
process (e.g. magnetic reconnection) taking place in a small 
region of the inner disk. This heating process has its own life- 
cycle that can be described with a rise time and a decay time 
of each component. 

The timescales for these processes are determined by the 
synchrotron cooling time Tjyn, which depends on the mag- 
netic field strength B and the wavelength A of the observed 
sy nchrotron radiation. F or the case of Sgr A* it is estimated 
as jGiUessen et al.ll2006l) : 

Tsyn w 8 (B/30G)"^''^ (A/2^m)'''Vin (9) 

The magnetic field strengths predicted by radiatively ineffi- 
cient accretion flow (RIAF) models for Sgr A* lie around 
20G, yi elding cooling times that are les s than the duration of a 
flare (Ouataert 2003; Yuan et al."2004'). Computational mod- 
els in general relativistic magnetohydrodynamics (GRMHD) 
suggest that in the close vicinity of a Kerr black hole the mag- 



netic field strength within the accretion disk decreases with 
radial distance. In contrast, higher spin-parameters lead to 
stronger magnetic fields (De Villiers et al. 2003). These de- 
pendencies are fairly complicated and presumably cannot be 
derived analytically. Since a full GRMHD simulation exceeds 
the scope of this work, we empirically determine the heat- 
ing and cooling timescales from the observed light curves of 
Sgr A*. 

It must also be mentioned that up to now, self-consistent 
GRMHD simulations can only produce hot spots with very 
short life ti mes, amounting to frac tions of the orbital period 
(r w 0.3P, Schnittm an et al.ir2006l) . Hence, there have been 
attempts to explain th e obser ved light curves with multi-spot 
models dEckart et al] l2008ah . Several hot spots distributed 
radially around the black hole can produce similar quasi- 
periodic signals in the light curves as a single source. How- 
ever, with the number of spots the flux amplification in the 
light curves and the deflection of the emission centroid is in- 
creasingly reduced. Unknown mechanisms that may stabilize 
hot spots in the accretion flow may cause longer life times, 
though. 

We implement the arc by adding up all the sphere's ray 
traced images of one complete orbit. This results in a lensed 
ring around the black hole, including its multiple images. As- 
signing a spatial and temporal brightness distribution to this 
ring allows us to simulate any desired hot spot. We use sim- 
ple Gaussian functions to describe the azimuthal brightness 
distribution of the arc and the temporal evolution of the entire 
hot spot. We leave the radial brightness distribution of the hot 
spot unaltered, i.e. it remains confined to a radial region of 
0.5rs diameter (see Fig. |7J- 

The spectral indices of sphere and arc remain to be speci- 
fied. Observations of the G C-flares suggest a con nection be- 
tween brightness and color (Gilless en et al.ll2006l) . When the 
flare is dim it shows a red spectrum, whereas in bright states 
it is rather white. Since in our model the sphere is brighter 
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Fig. 7. — A composition of ray traced images of a spliere on the ISCO 
about a Scliwarzscliild black hole in 70° inclination. The resulting hot spot 
is composed of a bright compact sphere in its center and an extended arc that 
follows a Gaussian brightness profile. Due to relativistic beaming, this profile 
appears modified. 

than the arc, motivated by the observations, we use a'^'^^"'' = 
and g""^ = -3, thus covering the range of observed values 
dEisenhauer et al. 2005a ; Gillessen et al..,2006; .Krabbe et al.i 
l2006HHornste'inet al.il2007l) ! 

3.3. Fitting observational data 

In this section we fit the parameters of our hot spot model 
to observational data from Sgr A*. To this end a x^' 
minimization is applied to fit the simulated light curves to the 
observed ones. However, for the fit of a flare light curve with a 
given periodicity we have to restrict ourselves to those model 
configurations with nearly the same orbital timescale. Thus, 
only grid-configurations with fixed r and a can be considered 
and leave us with four different inclinations per fit (our grid 
does not contain configurations with different values for r and 
a for a given orbital period). However, it is not the objective 
of this paper to estimate the spin-para meter of Sgr A* (this 
has akeady been done by 'Meyer et al.l f2006a,b)). We rather 
want to focus on the astrometric study of the hot spot model. 
Hence, for the parameters describing the hot spot itself, we 
can calculate optimized values. 

We used photometric measurements of two bright GC- 
flares, namely the flare from April 4 and from July 22, 2007. 
A detailed discussion of the fitting procedure is given in the 
following subsections. 

3.3.1. L-band flare from April 4, 2007 

Photometry — The light curve of this flare has three major, 
fairly broad peaks with quite a substructure (see Fig. [8]l. 
Additionally, four smaller peaks are found on the wings of 
the light curve. The power spectrum indicates a periodic 
timescale of about 23 min. Any details about the observa- 
ti on of this flare and the used instru mentation can be found 
in iDodds-Eden et al.l din preparationh . For the fit we choose 
grid-configurations with r = l.Srisco and a = 0.7, since they 
yield the closest orbital timescale of P = 22.6 min. 

In Fig. [8] the best fits for four different inclinations are 
shown. Each light curve constitutes seven orbital motions (ac- 
cording to the number of peaks) of a compact sphere (a^f*"™ = 



0) with an extended arc (a^"^ = -3). For the parameters of the 
hot spot we find the values given in Table [T] The decay time 
of the hot spot is slightly larger than its rise time, because the 
light curve is more flattened on the decaying flank. The arc 
provides quite a strong contribution to the total flux: The sur- 
face brightness at its center is only four times less than the sur- 
face brightness of the sphere. It is spread along the orbit, pos- 
sessing a full width at half maximum (FWHM « 2.355cra,-c) 
of more than one complete orbit. This value is required by the 
broad peaks in this flare fight curve. 



Parameters: 


Sphere 


Arc 


Rise time a,isc 


I.IP 


I.IP 


Decay time trjecay 


1.8P 


1.8P 


Brightness ratio 


4 


/ 1 


Extension (FWHM) 




1.41C 



TABLE 1 

Best wt parameters for the L-band flare from April 4, 2007. 

P IS THE orbital PERIOD AND C = Inr THE CIRCUMFERENCE OF THE 
HOT SPOT ORBIT. 



The reduced -values favor lower inclinations. However, 
one has to be careful with drawing conclusions from this for 
the following reason: The arc-width (Tare and the inclination / 
are, to a certain extent, degenerate parameters that both deter- 
mine the width and the height of the light curve peaks. Nev- 
ertheless, the implementation of an extended arc is necessary 
in order to fit the broad light curve peaks in the flare. The 
inclination only moderatly influences the flare shape, by far 
weaker than the arc- width does. 

The fit worsens towards the wings of the flare. By check- 
ing the spacing between the major peaks one finds that the 
outer peaks do not match the 23 min timescale as well as the 
inner ones. One could argue for either an acceleration or a 
deceleration of the orbital period, corresponding to shrinking 
or growing orbits, respectively. But, a close look at the spac- 
ing between the outer peaks does not support this reasoning, 
since their occurrence deviates from the orbital periodicity in 
both ways. 

Higher inclinations produce sharper sub-peaks on top of the 
major peaks. These are generated by the compact sphere and 
reproduce some substructure in the observed flare quite well. 
The Einstein ring peaks in the 90° -case are still compatible 
with the data, however, the sampling of the observed light 
curve is too sparse to constrain them. An unambiguous detec- 
tion of such a peak could directly suggest the occurrence of 
an Einstein ring, a high-order general relativistic effect. 

Having found a reasonable fit to the observed flare light 
curve it is of interest to extract its periodic properties by com- 
puting a power spectrum and comparing it to the one obtained 
from the experimental data. To thi s end we use a Lomb nor- 
malized periodogram (lLombllI976t IScargleiri982l) which is a 
standard tool for analyzing unevenly sampled time series. It 
is normalized by the variance of the overall flare light curve 
flux. 

The power spectra for the four best fit light curves from 
above are shown in the right inset of Fig. [S] Obviously, the 
23 min feature is nicely reproduced in the simulated data. 
Higher inclinations show more power in the periodogram due 
to the fact that the peaks in the corresponding light curves 
are more pronounced. The rise towards longer timescales is 
caused by the overall structure of the flare light curves, yield- 
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Fig. 8. — Left: Best fit light cui-ves for the L-band flare from Apiil 4, 2007 for four different inclinations (i = 20° (dotted), 50° (short-dashed), 70° (long- 
dashed), 90° (solid)). The gray error bars represent the measurement, at the bottom the residuals of the fits are drawn. Right: Power spectrum of these light 
curves (measurement in gray). 




ing a periodic timescale of the whole flare duration. Higher 
frequencies, though, are suppressed in comparison to the ob- 
served data. This is due to the fact that the source as well as 
the observing instrument generate intrinsic noise. While the 
power spectral density (PSD) of the instrumental noise is usu- 
ally Gaussian (i.e. constant), the one generated by the source 
may have a more complicated profile (like a broken power 
law, for instance). In this case the overall power spectrum of 
the observed flare shows a slight rise to longer timescales. 

In order to include the observational process into our sim- 
ulations we artificially add random noise to our data. We can 
also mimic the flux-integration at the detector of the telescope 
by averaging over several successive data points. Fig.|9]shows 
the light curves from above, treated with the said procedure. 
First, the light curves were integrated with an integration time 
of 45 seconds, which is the average value from the observa- 
tion. Then, a normally (Gaussian) distributed noise with a 
mean of zero and a standard deviation of 5 percent of the flare- 
peak-flux was superimposed. This value can be obtained by 
comparison of the mean error in the flare-data to its maximum 
flux-value. 

The corresponding power spectra are presented in the right 
panel of Fig. |9] The power in the high frequencies is raised 
and has a flat shape, pretty similar to the observed data (the 



fact that there is more power in the simulated periodograms 
at short timescales may be due to an overestimation of the 
observational errors). The peaks associated with the 23 min 
timescale now match up very well. The fact that we add Gaus- 
sian noise to the simulated data makes the power in this char- 
acteristic frequency decline. Thus, the observed light curves 
are consistent with the assumption of a clearly determinis- 
tic physical process (orbital motion) plus a superposed instru- 
mental Gaussian noise. 

Astrometry — Our code additionally provides simulations for 
the centroid motion of the flare. This is a direct prediction 
of the hot spot model and may be checked by future observa- 
tions. In Fig.[TO]the centroid tracks of the four best fit models 
from above are displayed. They seem to show strong rela- 
tivistic effects, however the axis labelling reveals the scaling. 
Due to the rather strong contribution of the extended arc, the 
centroid of emission does not move far, yielding deflections 
of about 1 .5 to 6 fias during one orbit. 

The centroid track is dragged to the approaching side be- 
cause of relativistic beaming; lensing effects shift it upwards. 
The higher the inclination the stronger are these effects, ex- 
cept for the 90°-case where symmetry avoids a y-direction 
shift. Due to the delay between primary and secondary im- 
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Fig. 10. — Centroid tracks of the best fit model for the L-band flare from April 4, 2007 for four different incHnations. The brightness evolution of the flare is 
indicated by the line thickness (the thicker the line the brighter the flare). A cross marks the beginning, an arrow the end of the flare. 



age, the centroid track slowly moves upwards, following the 
emission centroid of the primary image. A longer flare dura- 
tion would cause the centroid track to move down again, one 
can see that at the end of the simulation the motion is just at 
its turning point. 

Although it may be impossible to astrometrically resolve a 
flare like this in the near future, it is noteworthy that the smear- 
ing of the light source relatively strengthens the relativistic 
features in the centroid tracks. The motion of the overall cen- 
troid track is an indication of the secondary image, however, 
the shift is marginal. 

3.3.2. L-band flare from July 22, 2007 

Photometry — Another fairly bright flare from Sgr A* was 
caught on July 22, 2007 at the VLT under good observing 
conditions. There were a few short interruptions because of 
sky observations, causing some small gaps in the light curve 
data. The details about this observation can be found in 
iHaubois et al.l din preparation!) . Qualitatively, this light curve 
appears quite different from that of April 4. The major peaks 
are more prominent and the overall rise and decay of the 
whole flare light curve is suppressed. In particular, a wider 
separation of the peaks is noticeable, amounting to roughly 45 
minutes. In order to fit this flare one needs to assume a larger 
orbital radius of the orbiting material and a higher contribu- 
tion of a compact emission region. To this end we use a sim- 
ilar hot spot model as before and find the best fit parameters, 
reported in Table |2] We used configurations with r = 2.0risco 
and a = 0.52, yielding an orbital period of 45.4 min and sim- 
ulated 5 orbits. 



Parameters: 


Sphere 


Arc 


Rise time a-jise 


0.5P 


4.0P 


Decay time cr^^^.^y 


l.OP 


7.0P 


Brightness ratio 


86 


/ 1 


Extension (FWHM) 




1.24C 



TABLE 2 

Best fit parameters for the L-band flare from July 22, 2007. 
p is the orbital period and c = iwr the circumference of the 
hot spot orbit. 



The resulting best fit light curves are superimposed on the 
flare-plot in Fig.[TT] Because the sphere is dominant in this fit, 
the 90°-case shows very pronounced Einstein ring peaks. Un- 
fortunately, the strongest peak is not covered with data. The 
arc only makes a small contribution to the light curve flux and 
has the effect of broadening the wings of the major peaks. 
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Fig. 11 .— Best fit light curves for the L-band flare from July 22, 2007 for 
four different inclinations (/ = 20° (dotted), 50° (short-dashed), 70° (long- 
dashed), 90° (solid)). The gray eiTor bars represent the measurement, at the 
bottom the residuals of the fits are drawn. 

The arc's rise and decay time is so long that it barely varies 
in brightness. In contrast, the sphere must show stronger vari- 
ations in order to reproduce the observed peaks. Although 
its rise and decay times seem to be very short in terms of 
the orbital period, they are comparable to the corresponding 
timescales of the April flare. This concordance indicates a 
fundamental heating- and cooling-process of the flares which 
does not change from one event to the other (such as a mag- 
netic reconnection). 

A realistic simulated observation is shown in Fig.fT2lusing 
an integration time of 50 seconds and a Gaussian noise with a 
standard deviation of 6 percent of the flare-peak-flux (average 
values from the data). The corresponding power spectrum is 
also given in that figure. There is quite a good agreement be- 
tween the simulations and the observation. The power in the 
characteristic timescale (orbital period) exceeds the one in the 
April flare. This is due to the concentrated emission of the 
compact sphere in the center of the emission region, which 
produces more accentuated peaks in the light curves. Here, 
also the intermediate timescales yield a good agreement. The 
slightly too high power in the upper frequencies again sug- 
gests an overestimation of the errors in the measurement. 

Astrometry — The astrometric motion of the centroid, as it 
results from this fit, yields better prospects regarding future 
measurements than the April flare. As can be seen in Fig. [13] 
diameters of up to 40 /ias are reached in the centroid tracks. 
This is due to the bright compact component in this model and 
the larger orbital radius of the hot spot. Because the sphere's 
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Fig. 13. — Centroid tracks of the best fit model for the L-band flare from July 22, 2007 for four different inclinations. The brightness evolution of the flare is 
indicated by the line thickness (the thicker the line the brighter the flare). A cross marks the beginning, an arrow the end of the flare. 



brightness varies strongly, the centroid deflection grows first 
and then shrinks again, until it finally reaches its starting 
point. 

As for the light curves, observations of the centroid motion 
as they might be possible in future experiments can be simu- 
lated. For these we again integrate every 50 seconds and add 
Gaussian noise to both the x- and y-component of the cen- 
troid. Its standard deviation is chosen to be 10 /las, the de- 
sired resolution of ground-based interferometric instruments 
like GRAVITY. 

Fig. [14] displays the data points obtained from this proce- 
dure. In order to shrink the error bars, it is useful to average 
over several successive data points (here about 10). In the 
cases of 50° and 70° inclination the plots additionally contain 
the x{t) versus 0(f) motion of the centroid (shifted to x(0) = 0) 
to visualize its deflection more clearly. 

It appears that the centroid tracks at higher inclinations are 
harder to catch. This is due to the increasing beaming effects, 
dragging the centroid towards the approaching node. Nev- 
ertheless, significant deflections of the centroid can be mea- 
sured in most cases. A comparison of many different flare 
events could further constrain the configuration of hot spot 
and black hole. 

4. DISCUSSION 

The hot spot model we applied in this study is in good 
agreement with the observed light curves and periodograms 
from Sgr A* and bears good prospects for the successful de- 
tection of a centroid motion with the GRAVITY-experiment. 



However, some assumptions of this model still need to be jus- 
tified. 

Compared to ot her authors (iBroderick & Loebl l2006al : 
iMeyer et all 12006 b') we used a rather small and distinct hot 
spot component for the c entral emission region (following the 
estimate of r < 0.3rs by lGillessen et alJ (l2006h V The result- 
ing light curves and centroid tracks with o nly this component 
compare well to the results obtained by IBroderick & Loea 
(j2006a.b), except for the fact that the features are more dis- 
tinct due to the smaller spot size. 

The implementation of an arc is motivated by the broad sub- 
flare structures, which cannot be fitted well with a single spot, 
but an extended emission region. Thus, only the combination 
of these two components allows a reasonable fit to the ob- 
served data. We did not implement any disk component, since 
it does not contribute any dynamical information to the sim- 
ulated d ata and can be subtracted i n the optically thin regime 
(see e.g. lBroderick & Loebll2006bl) . 

Although relativistic hot spot models and positional dis- 
placements in the context of Sgr A* have been studied deeply 
in the recent past, observational data has never been applied to 
constrain the astrometry. In this respect, this work is unique. 

Furthermore it needs to be mentioned that the VLTl mea- 
surements will presumably be carried out in K-band. Since 
this study only covers L-band observations it might lead to 
a different picture. However, observations bv iHornstein et alj 
(.2007.) have shown that light curves from Sgr A* are well cor- 
related at both wavelengths. Also theoretical SED-models of 
the flares suggest a strong correlation between K- and L-band 
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Fig. 14. — Centroid tracks from Fig. 1131 superimposed by simulated data points as expected to be obtained with GRAVITY. The cases of 50° and 70° 
inclination additionally contain the x{t) versus </)(() motion of the centroid shifted to a:(0) = 0. 



flux due to the rather broad synchrotron bump created by a 
population of hea ted electrons in the magnetic reconnection 
(I Yuan et ani2003h . Therefore, the predictions for the astro- 
metric deflections from the light curves are independent of 
the wavelength band (K or L). 

It is being argued in recent publications which kind of phys- 
ical processes could generate the observed profiles in the PSD 
and how to distinguish betw een a pure red noise process and 
a distinct Keplerian signal. iDo et alJ (|2008) investigate the 
NIR flare light curves observed with the Keck telescope using 
Monte Carlo simulations. They claim to find no significant 
periodicity at any timescale probed and that periodic oscil- 
lations with amplitudes greater than 20 percent of the flare 



peak flux do not exist. Instead, their observations appear to be 
consistent with a red noise process with a power l aw slope be- 
tween 2.0 and 3.0. Similar results are attained by' Meyer et alJ 
(120081) who use data taken with both the Keck and the VLT, 
whereas feckart et akl ('2008b) interpret the variations of sub- 
structures in their recently observed flare data as variation of 
the hot spot structure due to differential rotation. 

The quality of the latest flare data is still not considered 
good enough to unambiguously rule out a random noise pro- 
cess in favor of a rotating hot spot model or vice versa. Never- 
theless, it is remarkable how well our simulated power spec- 
tra fit the measured ones. The basic assumption of a peri- 
odically moving object in our model does not produce a sig- 
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nificantly higher power in the orbital timescale than the ob- 
served characteristic timescale. Also the rise towards longer 
timescales below the orbital period is present. The reason for 
this is presumably the slightly non-sinusoidal nature of the 
flux-variations, as apparent from Fig.|2] 

Although our model is in a way idealized and more complex 
scenarios such as inwards spiraling spots and hydrodynamic 
instabilities in the accretion flow might also exist, there may 
still be a fair chance that single bright hot spots can dominate 
the emission from Sgr A* for several hours. In this case the 
hot spot model could be tested with astrometric observations 
as planned by the GRAVITY-experiment. 

5. CONCLUSIONS 

We have used Sgr A* flare observations in the L-band to test 
the hot spot model as a physical origin for the observed quasi- 
periodic oscillations. In order to predict astrometric properties 
of the flares we developed a model for the hot spot emission 
and applied general relativistic ray tracing simulations to vi- 
sualize its appearance for a remote observer. 

In the case of a confined spherical emission-region (small, 
constant source) various relativistic effects can be studied in 
detail. Both light curves and centroid tracks show strong devi- 
ations from the expectations of Newtonian physics, especially 
if the hot spot is close to the black hole and its orbit is highly 
inclined. In such configurations the occurrence of Einstein 
rings and multiple images can be inferred from substructure 
in the data. Furthermore, we presented a new approach to 
estimate the parameters of the hot spot model using astromet- 
ric information. However, this model is rather suited to study 
the fundamental properties of these parameters in a qualitative 
manner, not to fit the observations quantitatively. 

The use of polarizing filters for the analysis of centroid 
tracks in principle enhances the signatures caused by multi- 
ple images. However, due to the low degree of polarization in 
the flares and the little resulting flux in the individual polariz- 
ing channels, a measurement of this type is probably not very 



promising yet. 

A more advanced model for the hot spot assumes an elonga- 
tion of the compact source with a certain temporal evolution 
(extended, variable source). The parameters of this model can 
be adjusted in such a manner that the observed light curves 
from Sgr A* can be reproduced fairly well. An analysis of 
the simulated light curves in frequency space yields a good 
agreement with the observed data. 

In view of the observational data at hand (light curves) our 
model can only constrain the orbital period and the heating 
and cooling timescales unambiguously. Other parameters of 
the hot spot model, such as the orbital radius, the inclination, 
the arc-extension, or the spin of the black hole, are degenerate. 
Some of these degeneracies can be broken with astrometric 
data, such as the degeneracy in orbital radius, inclination and 
black hole spin. Nevertheless, the impact of the arc-extension 
and the brightness-ratio between sphere and arc remains to be 
disentangled from the orbital parameters. Future high quality 
polarimetric measurements could help in this respect. 

The centroid tracks from the best fit models are encourag- 
ing for future high resolution measurements. In particular the 
July 2007 flare announces the possibility to detect a centroid 
motion. The future interferometer GRAVITY at the VLTI 
will most likely be able to resolve proper motion of individual 
flares. Such a measurement would not only further constrain 
the hot spot model, but also probe spacetime in the strong 
field limit of general relativity. Thus, our simulations provide 
theoretical predictions from fundamental physics that can be 
tested in the near future. 
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